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Abstract 

We present a 'super symmetric' modification of the d-dimensional quantum 

rotor model whose ground state is exactly soluble. The model undergoes 

a vortex-binding transition from insulator to metal as the rotor coupling is 

varied. The Hamiltonian contains three-site terms which are relevant: they 

change the universality class of the transition from that of the (d+ 1) — to the 

d-dimensional classical XY model. The metallic phase has algebraic ODLRO 

but the superfluid density is identically zero. Variational wave functions for 

single-particle and collective excitations are presented. 
PACS: 05.30.-d,05.30.Jp,67.40.db,67.40.Vs 
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This paper discusses an exactly soluble modified version of the quantum XY model 

«o-[/E^- J£cos(^ - 9 j+s ), (1) 

3 3 j5 

where j is summed on sites of a hypercubic lattice and 5 is summed on near-neighbor 
vectors! This model is often used to describe the superconductor-insulator transition in 
granular superconductors and Josephson junction arrays.i The coupling constant J repre- 
sents the strength of the Josephson coupling between the order parameter phases 9i and 
9j on neighboring grains. The parameter U represents the charging energy of the grains. 
The boson (Cooper pair) number operator conjugate to the phase is the angular momentum 
fij = —idg.. While this is correctly quantized in integer values, it can be negative. Thus 
the model implicitly assumes a large background number uq of bosons per lattice site so 
that fij represents local deviations (positive or negative) from this mean (integer) value. We 
can view the cosine term as a mutual torque which transfers quanta of angular momentum 
(bosons) from one site to the next. Thus the quantum XY model is essentially equivalent 
to (i.e., in the same universality class as) the boson Hubbard model. For large U/J the 
ground state is a Mott-Hubbard insulator and for small U / J it is superfluid which exhibits 
off-diagonal long range order (ODLRO) in the phase field correlations (at zero temperature) 

Gij = (e i6i e-^) (2) 

for dimension d > 1 (and algebraic ODLRO for d — 1). The transition between the superfluid 
and Mott-Hubbard insulating states is continuous and is in the universality class of the d+1- 
dimensional classical XY modelJUil The extra dimension arises from the fact that in the path 
integral representation of the partition function, the Euclidean time interval < r < Tif3 
diverges at zero temperature. 

While the physics of this model is now completely understood, it has resisted exact solu- 
tion in all dimensions. It is interesting to consider a Jastrow-like variational wave function 



^o(0i, ■■■,9n) =exp 



.(hv(e l ,...,e N ) 



(3) 
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where A is a variational parameter and 

K = _j^ cos ^._^ j) (4) 

is the potential energy from eq. ([!])■ This form is motivated by the harmonic oscillator 

r? 1 

H = -^— + -Kx 2 (5) 

2m 2 V ' 



for which the exact ground state is 



ip(x) = exp 



-huo{^Kx 2 ) 



(6) 



The variational state in eq.(||) is thus in the spirit of the harmonic spin- wave approximation 
in which one expands the cosine term to second order in deviations from the classical ground 
state. The wave function is much better than this however because it obeys the correct 
periodicity under 9j — > 8j + 2n. This feature is crucial to the existence of (quantum) 
vortices in the ground state and hence allows for the possibility of a phase transition to the 
insulating state-physics which is completely missing from the spin-wave approximation. 

The purpose of the present paper is to examine the variational wave function of eq. (|3|) and 
to consider a Hamiltonian for which ip happens to be the exact ground state.i The question 
of what Hamiltonians have Jastrow wave functions for ground states has a long historyi'i 
and this question has recently been reexamined from a modern perspective by Kane, et al.i 
Except for special casesjll the generic requirement is that the Hamiltonian have three-body 
interactions of a particular form. We will see shortly that the analog of this for the present 
problem is three-site interactions. Kane, et al. argue that these three-body interactions 
are (perturbatively) irrelevant in the renormalization group sense. That is, they have no 
effect on the long-distance properties of the system (other than a trivial renormalization of 
the speed of sound). This result is perturbative because it neglects vortex-like excitations.^ 
While, strictly speaking, the above statements are true, they can be quite misleading if the 
ground state of the system undergoes a phase transition. We show explicitly below that 
the universality class of the transition is completely altered by the inclusion of particular 



three-site interactions. This is a non-perturbative effect precisely due to the role played by 
vortices in the (zero temperature) transition. The discussion below in readily generalized to 
any dimension, but for definiteness we consider only the case d = 2. 
We construct the desired Hamiltonian by defining the operators 

Qj = + 7u ? [sin {6j ~~ 6j+s) + sin {6] ~ 6j - s)] (7a) 

Qj = "V^ J- + £[rin (9 j - 9 j+s ) + sin (9, - 9^)}. (7b) 
It is readily verified using Eq. (|) that 

QAo = (8) 

for every j (we take A = 1 hereafter). 
The 'supersymmetric'ili Hamiltonian 

# = E Q)Qj (9) 
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is clearly positive semidefinite and therefore i/jq is an exact, zero-energy ground state of H. 



Using eq.(7a) we can write H in the form 



H = -uY,$ i - JE cos (%-^) 

+ ^ E Sin ~ Oi+S) Sin - e i+*)- 
U jSS' 

The first two terms on the righthand side are equivalent to the usual quantum XY model of 
eq.(|l]). The remaining term is a perturbation consisting of two- and three-site interactions. 
These terms represent the simultaneous hopping of a pair of bosons. The 5 = 5' terms give 
rise to a cos 28 coupling, a form which has been studied by Lee and Grinstein.0 

Now that we have the Hamiltonian and the exact ground state, let us examine the nature 
of the ground state as a function of the quantum fluctuation parameter U / J to see if the 
system undergoes a phase transition. The (unnormalized) probability distribution of the 
phase angles is 



p[9i, e 2 , ■ ■ ■ , e N ] = |^ol 2 = ex P | ^ cos (*i - Oi+s) } ( n ) 

which is identical to the Boltzmann factor for the classical 2D XY model with U/2J playing 
the role of dimensionless temperature. This model undergoes a Kosterlitz-Thouless phase 
transition at a critical temperature^ U/2J * 0.9. This is clearly a different universality 
class from that of the usual 2D quantum XY model which is known to be in the universality 
class of the 3D XY model. Thus we see that the three site terms are relevant to the transition 
(at least when they have the particular strength given in eq.(|T0"D). 

If we knew the exact ground state ■ ■ • , 6n) of the usual quantum rotor problem, 
then |$| 2 would of course define a 2D classical statistical mechanics problem. However, the 
fake classical Hamiltonian would necessarily contain long-range forces (in order to give the 
3D XY universality class in a 2D model!) . 

We are used to the notion that thermal fluctuations produce vortices. Here we see a nice 
illustration of the fact that even at zero temperature, vortices can be produced by quantum 
fluctuations. For U/2J > X£ T the largest amplitude configurations in the ground state 
contain free vortices and the spin-spin correlation function decays exponentially 

= (e ie 'e- i0 >) ~ e^'-^ 1 ^, (12) 

whereas for U/2J < T£ T the correlations decay only algebraically 

Gtf~ {Ri-Rjl^ (13) 

because vortices are confined. That is, virtual vortex-antivortex 'vacuum fluctuations' ap- 
pear but do not proliferate. 

Knowing the ground state exactly, we now find approximate excited states of the Hamil- 
tonian. Let us assume an excited state loosely analogous to the Feynman-Bijl formS 

* ? =p**o- (14) 

As discussed below, the form of p$ depends on the choice of excitation. If the excited state 
is orthogonal to the ground state, i.e., 



(tf,H*o> = (15) 

then we can use the variational principle to find an upper bound on the excitation energy 
(taking advantage of the fact that the ground state energy vanishes) 

where (^fg\H\^f 9 ) = /(g) is the 'oscillator strength', and (^ q \^/ q ) = s(q) is the 'static 
structure factor'. Writing is terms of ty : 

f(q) = (9dH\9 9 ) 

= (* |p-^pd*o), (17) 

and writing the Hamiltonian in terms of Q's and using the fact that Qj^o — 0, 

f(q) = (^o\[p- < !,Q]][Q 1 ,pi\^) (18) 

By substituting an explicit form of in the above results, f(q) can be calculated imme- 
diately. 

First we consider an excited state wave function which describes a single-particle excita- 
tion: 

*?=&J*o (19) 

where: 



bi=^Ye i * R ie i0 * (20) 

q Vn^ 



i 

j ■ . 

is the spatial Fourier transform of the operator that adds a unit of angular momentum at 
site j. Using Eq. (|20"D 



[QiM = \f u / Ne n ^ j (21) 

it follows that = {9o\b-cfHb^o) = U. The static structure factor is given by the spin 
susceptibility of the classical XY model at wave vector q 



s (?) = jrJ2 ei€ ' {Mi ~ nj) (^o\e i{ei ' 6j) \%) (22) 

ij 

which we know from Kosterlitz-Thouless theory: 

\Rj — Rj\ ^ T < T r 
(e^ e -^) ~ " (23) 

where T c is the critical temperature, or critical coupling (U f J) c in our case, with 77 ranging 
from to 1/4 as temperature varies from to T c . In the spin- wave approximation r\ = -^j 
for our model. Substituting Eq. (p3|) in Eq. ([22]), and changing summations to integrals for 
an infinitely large system with U/J below the 
critical point 

s{q) / r^ +1 J (qr)dr ~ g" 2 ^, (24) 

2n Jo 

where r = \Ri — Rj\ and since f(q) = U 

A(q) ~ iL ~ itf"". (25) 
For U / J above the critical point and g£ <C 1 

s(g) = — / drre- r/ tj (qr) ~ £ 2 (26) 

27T JO 

and hence the quantum system has an excitation gap (within the single mode approximation) 

A(q = 0)~UC 2 . (27) 

Using the Kosterlitz-Thouless theoryE! prediction for the correlation length we have 

A(0)~exp{-2^-(i)j- 1/2 } (28) 

where b is a positive constant. 

Figs. (|l]-0) illustrate the basic features of the excitation energy. In the thermodynamic 
limit the system is gapless below the critical point and has a gap which rises from zero above 
the critical point with the essential singularity characteristic of the KT transition. In the 



ordinary quantum XY model, the Bogoliubov process mixes the single-particle and density 
excitations to produce a linearly dispersing collective Goldstone mode in contrast to the 
uj ~ q 2 dispersion of free bosons. We see here from Eq. ( |25|) a curious contrast to the generic 
behavior. The collective mode dispersion u ~ q 2 ~ v gradually stiffens with increasing U/J 
but never becomes linear since r\ < 1/4 below the transition. 

We note that right at the critical point the q = single-particle energy only vanishes in 
the thermodynamic limit L — > oo since the classical XY model susceptibility obeys s(q) ~ 
/ d 2 r r _1//4 ~ L 7 / 4 . Thus we expect L 7 / 4 A(0) to be scale-invariant (i.e. independent of L) 
at the critical point, provided L is large enough. Using data from lattices with 8 < L < 24 
we found the scale-invariant point to be U/J ~ 0.905. On general renormalization group 
grounds we expect logarithmic corrections to scaling for small L in the 2D XY model. We 
include these corrections in Fig. ^ where we plot L 7//4 A(0)/[1 + 1/(2 In L + 4.5)] vs U/J. We 
again find the critical value U/J ~ 0.905 but the scaling now works well all the way down 
to L — 4. Our value for the critical coupling is close to, but somewhat above, the value of 
U / J ~ 0.895 found by Olsson and Minnhag en0 using the scaling of the superfluid density. 

We turn now to a study of the collective density mode excited state by taking to be 
the Fourier transform of the number density 

P 9 -=E^ l H^) (29) 

i 

so that 

= Tf E ^ ( sin & - d i+s) + sin (°i ~ 0i-*)}*o- (30) 
The static structure factor is 

= E eVft-fy (* | [sin (6 r - 9 i+s ) + sin (6, - 6^)} 

ijSS' 

[sin (9j - 9 j+s ,) + sm (9j - %-a0]|*o)- (31) 

Following the steps in the calculation of oscillator strength for the single-particle model, we 
find f(q) for the density wave state to be 
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/(g) = (* |_ £{(1 - e^)(l - e-^') cos - i+5 ) cos - 

U jS5> 

+ (1 - e-^)(l - e~^') cos - Ojs) cos - i+ e, to ,) 
+ (1 - e l ^)(l - e**) cos (9 j - 6 3+5 ) cos (9, - O^it*') 

+ (1 - e-^)(l - e**) cos - 9j- s ) cos - 9^eita')}\^o). (32) 

We performed Monte Carlo simulations of the 2D XY model to find the excited 

state energy < f(q)/s(q) for systems of finite size. Though one expects some physical 
connection between single-particle and density-mode approximations, our results for the two 
models are quite different. We found out that unlike single-particle excitation, the density 
wave is nearly dispersionless. Consequently, the results do not change significantly with the 
system size in the latter case. The excited state energy vs. q is shown in Fig. (|j), with 
q in (1,0) direction. Fig. ([5]) is a plot of vs. U/ J for the smallest non-zero allowed vector 
on a 32 x 32 lattice, i.e., ^§f,0^. We notice that the excited state energy for this model 
increases gradually with coupling U/J, contrary to the single-particle case, where the energy 
is close to zero below the critical point, and then abruptly increases. 

We conclude with some comments on additional curious features of this model. The 
model is readily extended to include exact solutions for arbitrary random bond strengths 
Jj j( 5 and frustration vector potential A^s- One might imagine that since the ground state en- 
ergy is identically zero independent of the disorder realization, one could compute ensemble- 
averaged correlation functions without having to invoke the replica trick.0 This is not pos- 
sible however since the norm of the ground wave function 

does depend on the disorder. A second consequence of the ground state energy being zero 
for all Aj t s is that the superfluid density is identically zero at T = even though the system 
exhibits (algebraic) ODLRO below the critical value of U/J. 
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FIGURES 

FIG. 1. Excitation energy A^ as a function of q for U / J = 0.8 and 0.9. The ( x ) are for q in the 
(1,0) direction and the (<0>) are for q in the (1,1) direction. For small q, the single particle model 
is isotropic. 

FIG. 2. Single particle excitation energy vs. coupling at q = 0, for different system sizes. Notice 
that Ag- — > as the system size increases for U/J below the critical point. 

FIG. 3. Estimation of critical coupling for the single-particle model. L 7 / 4 A(0) is independent 
of system size L at the critical point according to the Kosterlitz-Thouless theory. The logarithmic 
factor is a correction to scaling (see text). The parameter C ~ 4.5 was adjusted to obtain the best 
scaling. 

FIG. 4. Excitation energy vs. q for single-mode density-wave approximation, q is in the 
(1,0) direction. The value of the coupling U/J is shown next to each curve. 

FIG. 5. The density wave excitation energy e^* increases gradually as a function of U/J in the 
single-mode approximation. The results are for q = ^|f,0^, the smallest non-zero allowed wave 
vector on a 32 x 32 lattice. 
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